On the ACDM Universe in f{R) gravity 
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Several different explicit reconstructions of f(R) gravity are obtained from the background FRW 
expansion history. It is shown that the only theory whose Lagrangian is a simple function of the 
Ricci scalar R, that admits an exact ACDM expansion history is standard General Relativity with 
a positive cosmological constant and the only way to obtain this behaviour of the scale factor for 
more general functions of R is to add additional degrees of freedom to the matter sector. 
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INTRODUCTION 



After more than one hundred years, General Relativity 
(GR) is still considered to be the best fundamental theory 
for the description of the gravitational action. When ap- 
plied to cosmology, assuming homogeneity and isotropy 
(encoded in the Friedmann-Laimatre-Robertson- Walker 
metric), together with a fluid description of baryons. Cold 
Dark Matter (CDM) and radiation, GR gives rise to a set 
of field equations which when solved produce the simplest 
expanding cosmology - the Friedmann model, governing 
the dynamics of the cosmological scale- factor a{t). This 
model has been remarkably successful, giving for example 
the correct light element abundances and explaining the 
origin of the Cosmic Microwave Background Radiation 
(CMBR). In the last two decades, however, advances in 
observational cosmology appears to suggest that if one 
wishes to retain the FLRW metric, the universe must 
have undergone two periods of accelerated expansion. 
The first period of acceleration, known as the inflation- 
ary epoch is needed to explain the flatness problem and 
the near-scale invariant spectrum of temperature fluctu- 
ations observed in the CMBR, while the second period 
explains the dimming of distant type la supernovae rel- 
ative to Einstein-de Sitter universe model. In order to 
explain these periods of acceleration, the strong energy 
condition (p -l- 3p > 0) needs to be violated. In the case 
of inflation, this is achieved by introducing a dynamical 
scalar field, while the present day acceleration is most 
easily explained with the introduction of a cosmological 
constant. The resulting description of the Universe, in 
which the recent expansion history is driven by a cos- 
mological constant and ordinary matter is dominated by 
a CDM component hasbecome known as the ACDM 
or Concordance Model [1]. Unfortunately, this beauti- 
fully simple phenomenological model, which appears to 



fit all currently available observations (Supernovae la [2| , 
CMBR anisotropics [^IjLarge Scale Structure formation 
[J] , baryon oscillations [5( and weak lensing [6|) is affected 
by significant fine-tuning problems related to the vacuum 
energy scale and therefore it is important to investigate 
alternatives to this description of the Universe. 

Currently, one of the most popular alternatives to the 
A CDM model is based on modifications of the standard 
Einstein-Hilbert action. This is due to the fact that these 
changes naturally admit a phase of late time accelerated 
expansion (an early Universe infiationary phase is also 
possible [7|). In this way Dark Energy can be thought 
of as have a geometrical origin, rather than be due to 
the vacuum energy or additional scalar fields which are 
added by hand to the energy momentum tensor. As a re- 
sult, the cosmology and astrophysics of modified gravity 
is currently an extremely active area of research (see the 
recent reviews [8| and references therein). 

One of the simplest extensions of GR is based on grav- 
itational actions which are non-linear in the Ricci cur- 
vature R and/or contain terms involving combinations 
of derivatives of R |9l4lq . An important feature of these 
theories is that the field equations can be written in a way 
which makes it easy to compare with GR. This is done 
by moving all the higher-order corrections to the curva- 
ture onto the RHS of the field equations and defining an 
"effective" source term, often described as the curvature 
fluid. Once this has been done, it is easy to see that 
how a change of sign in the deceleration parameter of 
the FLRW cosmology can occur, leading to a period of 
late-time acceleration. 

Studies of the physics of these theories is however ham- 
pered by the complexity of the field equations, making 
it difficult to obtain both exact and numerical solutions 
which can be compared with observations. Theses prob- 
lems can be reduced somewhat by using the theory of dy- 
namical systems 16^, which provides a relatively simple 



method for obtaining exact solutions and a (qualitative) 
description of the glob al dynamics of these models for a 
given f{R) theory [l3|. 

Another useful approach is to assume that the expan- 
sion history of the universe is known exactly, and to in- 
vert the field equations to deduce what class of f{R) the- 
ories give rise to this particular cosmological evolution. 

This has been done recently for exact power-law solu- 
tions for the scale factor, corresponding to phases of cos- 
mic evolution when the energy density is dominated by a 
perfect fluid. It was found that such expansion histories 
only exist for i?" gravity 18|. A more extensive analy- 
sis of reconstruction methods has been carried out in [l9| 
to obtain theories which give an approximate description 
of deceleration-acceleration transitions in cosmology and 
in |20| a powerful approach to reconstruction based on 
standard cosmic parameters instead of a time law for the 
scale factor was introduced. 

In this paper we perform a number of explicit re- 
constructions which lead to a number of interesting re- 
sults. We find, for example, that the only real valued 
Lagrangian f{R) that is able to mimic an exact ACDM 
expansion history for a universe filled with dust-like mat- 
ter is the Einstein-Hilbert Lagrangian with positive cos- 
mological constant. This does not mean that f{R) grav- 
ity is incompatible with an exact ACDM expansion his- 
tory. In fact we further show that in a universe filled 
with a both a minimally-coupled non-interacting mass- 
less scalar field and dust-like matter, a theory of gravity 
can be found which exactly mimics the ACDM expansion 
history, making it impossible to distinguish it from GR 
at the level of the FLRW background. Moreover, number 
of realistic models may also mimic late-time acceleration 
epoch approximately. 



FIELD EQUATIONS FOR HOMOGENEOUS AND 
ISOTROPIC SPACETIMES IN f{R) GRAVITY 

We consider the following action within the context 
of 4-dimensional homogeneous and isotropic spacetimes, 
i.e., the FLRW universes with negligible spatial curva- 
ture: 



A = 



1 



d^x. 



'^[f{R)+£,n] 



(1) 



where R is the Ricci scalar, / is general differentiable (at 
least C^) function of the Ricci scalar and £„i corresponds 
to the matter Lagrangian. Units are chosen so that c = 
SttG = 1. 

The field equations for homogeneous and isotropic 
spacetimes are the Raychaudhuri equation 



3H + 3H^ = - 



1 



P + 3p+f-f'R + 3Hf"R 



2/' 
+3f"'R^ + 3f"R , 



(2) 



where H is the Hubble parameter; the Friedmann equa- 
tion 



iH^ 



1 



p^M_l-3Hf"R 



(3) 



the trace equation 

3Rf" ^ p-3P + f'R-2f- 9Hf"R - 3f"'if ; (4) 
and the energy conservation equation for standard matter 

p = -e{p + p). (5) 

Combining the Friedmann and Raychaudhuri equations, 
we obtain 



R = 6H + UH^ , 



(6) 



which is the usual definition for the Ricci scalar for ho- 
mogeneous and isotropic flat FRW spacetimes. 

The Raychaudhury equation can be obtained by 
adding the Friedmann equation with it's time deriva- 
tive (using the energy conservation equation and the def- 
inition of the Ricci scalar). Therefore any solution of 
the Friedmann equation is automatically a solution to 
the Raychaudhuri equation and hence it is sufhcient to 
solve the Friedmann equation to reconstruct the theory 
of gravity. 



RECONSTRUCTION OF A f(R) THEORY THAT 
ADMITS AN EXACT ACDM MODEL 



As we know, the present day observations suggest that 
the variation of the Hubble parameter with the redshift 
is sufficiently well described by the relation 



i?W = \/f(l + -)^ + f 



(7) 



where po ^ is the matter density (which consists of the 
observed and the cold dark matter) and A is the cosmo- 
logical constant. In what follows we try and construct 
theories (belonging to the class of f(R) gravity), that ex- 
actly mimic the above expansion history. 

From equation ([T]) we see that the time derivative of 
the scale factor a{t) can be given as 



Po 
3a 



A 
"3 



(8) 



Here we have used the usual definition of the redshift 
1/a =1-1-0. From the above equation we can immedi- 
ately calculate the second derivative of the scale factor, 
which is given by 



1 



(«').a 



2Aa^ 



Po 



6a2 



(9) 



We know for a flat FLRW universe, the Ricci scalar is 
defined by 



-R-6( - + — 
a a'^ 



(10) 



Using equations ([U and dU in equation (fTO|) we obtain 
the Ricci scalar in terms of the scale factor as 



R{a) = 



4Aa3 + po 



(11) 



We would like to emphasize here that till now we haven't 
a priori assumed any specific theory of gravity. Equa- 
tion ([5]) is obtained by observations while equation ([TU]) 
is a purely geometrical result for flat homogeneous and 
isotropic spacetimes, independent of the theory of grav- 
ity. 

We can now invert equation (flUl) . to write the scale 
factor in terms of the Ricci scalar 



a{R) = 



Po 



i?-4A 



(1/3) 



(12) 



We note that, since the scale factor has to be real, we 
considered only the real root of equation (jlOp , discarding 
the other two complex conjugate roots. Also the positiv- 
ity of the scale factor implies that the Ricci scalar reaches 
the value 4A, asymptotically in an infinite time. From 
equation (|12p we can calculate the Hubble parameter and 
the time derivative of Ricci scalar in terms of the Ricci 
scalar and these are given by 



HiR) = 



1 



po 



A 



a{R)y 3a{R) 3 
R = RAaiR))^ 



pa 



3a{R) 



A 
3" 



(13) 

(14) 



Similarly, using the energy conservation equation we can 
write the matter density of the universe in terms of scale 
factor 'a', or alternatively p{a{Rj). 

To investigate which functions f{R) exactly mimic the 
ACDM expansion history, we substitute all the above 
quantities written as a function of the Ricci scalar into 
the Friedmann equation, obtaining 



3(i?-3A)(i?-4A)/"(i?) 
p{R) 



|-3AJ/'(i?) + i/(i?) 



0. (15) 



Since this equation has to be satisfied for all times (which 
imply for all R > 4A G 3ff), this becomes a differential 
equation for the function f{R) in i?-space. It is easy to 
see that the above equation is an exact inhomogeneous 
hypergeometric equation for the variable x = —3 + i?/A, 
with p{R) as the inhomogeneous term. The solution to 



the homogeneous part is given by 






(16) 



where a± = (-7± y73)/12, /3± = (-11 ± \/73)/12 and 
Ci,2 are arbitrary constants of integration. 

Let us now analyze the solution carefully. The two fi- 
nite poles of the hypergeometric equation are at i? = 3A 
and R = 4A respectively. Since the allowed range of R is 
i? > 4A, one of the pole is out of the range and the other 
is at the boundary. However we see that in this range 
X > 1. We know that the convergence of a hypergeomet- 
ric function for the variable 'x' is guaranteed if |a;| < 1, 
and otherwise the function is either divergent or complex 
valued. Indeed one can check explicitly that both the so- 
lutions of the homogeneous equation are complex valued 
for R > 4A. Hence to ensure a real valued function f{R), 
we must choose both the arbitrary constants Ci^2 to be 
zero. This is interesting as it shows that there cannot be 
any real valued function of Ricci scalar that can mimic a 
ACDM expansion for a vacuum universe. 

Therefore from now on we only consider the particular 
solution for the given inhomogeneous term p{R). Let us 
now explicitly reconstruct the theories for which a given 
matter field would obey a ACDM expansion history. We 
would like to note here that a different reconstruction 
approach related with local tests in modified gravity was 
considered in 1241. 



Reconstruction for dust-like matter 

Supposing the Universe is filled with dust-like matter 
{w = 0). From the energy conservation equation we have 



PW = 71 



p{R) = i? - 4A 



(17) 



Substituting in the Friedmann equation (jlSp . we get the 
particular solution as 



/(i?) - i? - 2A 



(18) 



which is the well known Lagrangian for general relativity 
with a cosmological constant. This result is interesting 
as it proves that the only real valued Lagrangian f{R), 
that can mimic an exact ACDM expansion history for 
a universe filled with dust-like matter, is the Einstein- 
Hilbert Lagrangian with positive cosmological constant. 
It is also important to note, however, that this is not 
the case, if we put A = 0. In that case the general solu- 
tion of the Friedmann equation is 



f{R)^R + C\R"+ +C2R"- 



(19) 



We see that the function is real valued for positive Ricci 
scalar and hence there exist classes of real valued function 
f{R), other than GR, that can mimic a dust- like expan- 
sion history without the cosmological constant. But even 
a very small value of the cosmological constant would 
break this degeneracy, and in that case the theory must 
be GR. 



Reconstruction for perfect fluid with equation of 
state p = — l/3p 

A perfect fluid with an equation of state p = —hp is 
physically interesting as it lies in the boundary of the set 
of matter fields that obey the strong energy condition. 
In GR such fluids give rise to a Milne Universe which 
is a coasting universe, and the Ricci scalar is propor- 
tional to the square of the Hubble parameter. However 
as we would see, even this kind of fluids can also mimic 
a ACDM Universe in higher order gravity. 

Using the equation of state in the energy conservation 
equation and supposing the present density of the fluid 
is pf we get 



p{a)^^^p{R)^[pf{R~4A)r/' 



(20) 



Now solving the Friedmann equation (jlSp , the particular 
solution is 



(21) 



f{R) = p{R - 4A)2/^ , 
where /U is a constant depending on pf. 

Reconstruction for multi-fluids 



Let us now consider that along with dust-like matter, 
a non-interacting stiff fluid is also present in the universe 
and their present densities are po and ps respectively. 
This scenario is also possible if we have non-interacting 
minimally coupled massless scalar field with the dust- 
like matter. In this case, from the conservation equation, 
total matter density is 

P(^) = ^ + ^ ^ Pi^) = (i?-4A) + ^(i?-4A)2 . (22) 
a-' a" Pq 

Substituting this into the Friedmann equation (|15p. we 
get the following particular solution: 



f{R)^fiiR + P2R^ + t^3, 



(23) 



where pnin = 1..3) are constants depending on po, p^ and 
A. We therefore conclude that if the universe is filled 
with minimally coupled non interacting massless scalar 
field with dust-like matter, then the theory of gravity 
described above would exactly mimic a ACDM expan- 
sion history and it is impossible to distinguish this from 
GR with present cosmological observations for the back- 
ground FLRW level. 



Reconstruction for non isentropic perfect fluids 

Non-isentropic perfect fluids are generally described by 
the equation of state 



P = Hp,a) 



(24) 



Using the energy conservation relation we get the re- 
quired differential equation for p as 



Pia).a 



-{p + h{p,a)) 



(25) 



In general the above equation may not admit a closed 
form solution. However the calculation gets much simpler 
if /i(p, a) is a separable function of the form 



/i(/9, a) — w{a)p 



(26) 



that is, the barotropic index of the fluid changes with 
time. In this case we can directly integrate equation (j25p 
to get 



p{a) = exp 



l + w(a) 



da 



(27) 



We then substitute this into the Friedmann equation to 
get the required theory of gravity. 

As a specific illustration, lets assume that the time 
dependent barotropic index is given by 



w{a) 



27 — va^ 
7 + va-^ 



(28) 



where 7 and v are constants. As we can see in early 
times the fluid has positive pressure and at late times 
this behaves like a cosmological constant. Using ([27|) 
and (|12|) , we get p in terms of the Ricci scalar as 



p[R) 



{vpo + jR- AAjf 
^1 



(29) 



Substituting this into the Friedmann equation (jlSp . we 
obtain the particular solution to be 



f{R) = PiR + P2R^ + P3R^ 



Pi 



(30) 



where Pn{n ~ 1..4) are constants depending on poj 7i ^^ 
and A. 

Let us consider another form of non-isentropic perfect 
fluids where the equation of state is given by 



p — wp + h{a) 



(31) 



This equation of state is physically interesting, due to the 
presence of the tuning function h(a), which can compen- 
sate the effect of the fourth order gravity at intermedi- 
ate times. In this way we can have a matter dominated 
Friedmann like epoch (necessary for structure formations 



in the Universe), followed by an accelerated expanding 
phase. 

Integrating the energy conservation equation (^5]) . we 
obtain 



p(a) 



- / 3a(2+3"')/j(a)da + Ci 



^-3(i+w) ^ ^32) 



where Ci is an arbitrary constant of integration. Thus, 
we see that the effective density is that of an isentropic 
perfect fluid together with a time dependent cosmological 
term. We can always tune this term to obtain a fit with 
cosmological observations. 

As a specific example let us consider h{a) = a^^^ and 
w ~ 0. This then implies that the matter field in the 
Universe is dust along with a time dependent cosmologi- 
cal term which diverges at the Big Bang singularity and 
goes to zero sufhciently fast at the later epochs. 

Integrating the energy conservation equation we obtain 

p{a) = ^ + -^ => p{R) = i?-4A+^(i?-4A)4 . (33) 
a-' a^^ po 

Solving the Friedmann equation, we obtain the following 
particular solutions 



f{R) = piR + ti2R^ + p:iR^ + piR'^ + M5 



(34) 



where pn {n — 1..5) are constants depending on p^, pi 
and A. 



RECONSTRUCTION OF APPROXIMATE ACDM 
MODELS 

Let us now consider a useful technique proposed in 
Ref.[2l|, where the Friedmann equations are written as 
functions of the number of e-foldings instead of the time, 
A^ = In — . In such a case, the Hubble parameter is given 
in terms of A^ as: 



H = g{N)=g[-\n{l + z)) . (35) 

Then, the first FLRW equation for a flat universe yields, 

cPfiR) 



= -^G {N [R)) (AG' (N {R)) + G" [N [R))) ■ 



3G{N{R)) + Ig'{N{R))]& 



dR^ 



fiR) 



E 



PiQf^O 



-3{1+Wi) 3(i+Wi)N{R) 



(36) 



where G{N) = g {N) = H^, and the Ricci scalar can 
be related with A^ via i? = 12G{N) + 5G'{N). Then the 
equation ([55)) can be resolved for a given Hubble expan- 
sion rate, and the corresponding f{R) theory is recon- 
structed. 



Cosmological solutions in f{R) gravity with the 
presence of an inhomogeneous EoS fluid 

We consider now a Universe governed by some spe- 
cific f{R) theory in the presence of a perfect fluid, whose 
equation of state is given by. 



p^w{a)p + Cia) 



(37) 



where w{a) and C(a) are functions of the scale factor 
a, which could correspond to the dynamical behavior of 
the fluid and to its viscosity. Let us write the FLRW 
equations for f(R) as following. 



2H + 3H^ = -{p'+Pf^j,)) 



(38) 
(39) 



where p' = jif^ and p' = ,,f^, . The pressure and en- 
ergy density with the subscript f{R) contains the terms 
corresponded to f{R) and are defined as 

PfW = jj {^^ 3i/i?/") (40) 

PHR) = J (^'/"' + 2HRr + Rf" + !(/ - i?/')) (41) 

Then, by combining both FLRW equations, and using 
the equation of state defined in ([37]), we can write 

C(a) = (w{a)pfi^R) - pf^n) -2H- 3(1 + wia))H^) f'{R{a)) 

(42) 
As C(a) just depends on the Hubble parameter and its 
derivatives, for some specific solutions, any kind of cos- 
mology can be reproduced. Let us consider the example. 



—H — GpC 



Gqa"" 



(43) 



where Gp and Gq are constants 



We can check that in 
this solution, the first term in the r.h.s. corresponds to 
a fluid with EoS Wp = —I + c/3 > — 1 , while the second 
term, it has an equation of state Wq — —I — d/3 < — 1, 
which corresponds to a phantom fluid. We could consider 
a viable f{R) function proposed in [2 25, which is given 

by 



FiR) 



aR 



m+l 



;3i?" 



1-f 7i?' 



(44) 



This function is known to pass the local gravity tests and 
could contribute to drive the Universe to an accelerated 
phase. Then, by introducing ([^5]) and (H^ in the ex- 
pression for ^(a) in ([42]) . we obtain the equation of state 
for the inhomogeneous fluid that, together with f{R), re- 
produces the solution ([43l) . which for c = 3 and d > 
reproduces the ACDM model, and probably drives the 
Universe evolution into a phantom phase in the near fu- 
ture. 



CONCLUSION 

In this paper we have extended a the reconstruction 
programme for f{R) gravity obtaining several interesting 
expHcit reconstructions. In particular we find that the 
only real valued Lagrangian f{R) that is able to mimic 
an exact ACDM expansion history for a universe filled 
with dust-like matter is the Einstein-Hilbert Lagrangian 
with positive cosmological constant. This does not mean 
that f{R) gravity is incompatible with an exact ACDM 
expansion history, only that one has to extend the the- 
ory somewhat in order that this possibility can be re- 
alized. For example in a universe filled with a both a 
minimally-coupled non-interacting massless scalar field 
and dust-like matter, a theory of gravity can be found 
which exactly mimics the ACDM expansion history, mak- 
ing it impossible to distinguish it from GR using measure- 
ments of the background cosmological parameters. It is 
then an interesting problem to probe, how the perturba- 
tions in these modified theories can break this degener- 
acy, by predicting different structure formations, growth 
factor or cosmological gravitational waves from GR, that 
can be experimentally verified [23|. In fact the mat- 
ter power spectrum for power law f(R) theories satisfies 
the requirement for scale invariance and has distinct fea- 
tures that can be detected by combining future cosmic 
microwave background (CMB) and large scale surveys 
(LSS) data. Moreover, it remains the number of realistic 
modified gravities which reproduce this late-time cosmic 
acceleration era approximately, where essentially any de- 
gree of accuracy may be achieved. In that case, it is then 
required to have local constraints, like the Solar System 
or Post Newtonian constraints (see for e.g [24| and the 
references therein, which rules out the Lagrangians which 
grow indefinitely as the Ricci scalar vanishes) to have a 
modified gravity theory that works for both local and 
cosmological scales. 

This work was supported by a MICINN (Spain)-NRF 
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